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Abstract: 

Sum rules for twist-3 distributions are reexamined. Integral relations between twist-3 and 
twist-2 parton distributions suggest the possibility for a 5-function at x — 0. We confirm 
and clarify this result by constructing h L and h\ (quark-gluon interaction dependent part 
of hi) explicitly from their moments for a one-loop dressed massive quark. The physics of 
these results is illustrated by calculating hi(x,Q 2 ) using light-front time-ordered pQCD to 
0(as) on a quark target. A 6~(x) term is also found in e(x,Q 2 ), but not in gT(x,Q 2 ), to 
this order in 0(as)- 



1 Introduction 



Ongoing experiments with polarized beams and/or targets conducted at RHIC, HERMES 
and COMPASS etc are providing us with important information on the spin distribution 
carried by quarks and gluons in the nucleon. They are also enabling us to extract informa- 
tion on the higher twist distributions which represent the effect of quark-gluon correlations. 
In particular, the twist-3 distributions gT(x,Q 2 ) and hi(x,Q 2 ) are unique in that they 
appear as a leading contribution in some spin asymmetries: For example, qt can be mea- 
sured in the transversely polarized lepton- nucleon deep inelastic scattering and appears 
in the longitudinal-transverse spin asymmetry in the polarized nucleon-nucleon Drell-Yan 
process ||. The purpose of this paper is to reexamine the validity of the sum rules for these 
twist-3 distributions. 

The complete set of the twist-3 quark distributions in our interest are given as the 
light-cone corelation functions in a parent hadron with momentum P, spin S and mass M: 

^e^(PS\mri^n)\ Q i\PS) 

= 2 [^(ar, Q V(S • n) + g T (x, Q 2 )S^ + M 2 g 3 (x, Q 2 )n»(S ■ n)] , (1.1) 
^e ax (PSmO)a^^(Xn)\ Q ,\PS) = 2 \h 1 {x,Q 2 )(S i y - S^)/M 

+h L (x, Q 2 )M{p^n u - p v n^)(S ■ n) + h 3 (x, Q 2 )M{S^n u - Sin")] , (1.2) 

^ x *(PS\m^n)\ Q i\PS) =2Me(x,Q 2 ), (1.3) 

where two light-like vectors p and n are introduced by the relation p 2 = n 2 = 0, n + = p~ = 0, 
PM = _|_ Mi n M an d 5/1 i s decomposed as = (S ■ n)p" + (S ■ p)n" + Sf_. The variable 
x represents the parton's light-cone momentum fraction and each function has a support 
for x on [—1,1]. The anti-quark distributions g ljT3 (x, Q 2 ), hi^^(x,Q 2 ) are obtained by 
the replacement of ip into its charge conjugation field Cip T in (|1.1|) - (|1.3|) and are related to 

the quark distributions as g^r^ix, Q 2 ) = gi,T,z{—X-, Q 2 ) and hi,T,3(%, Q 2 ) = — ^i,r,3( — x i Q 2 )- 
The sum rules in our interest are obtained by taking the first moment of the above relations. 
For example, from fll.lQ , one obtains 

(psmoh^m\ Q ,\ps) 

= 2 jf 1 dx [ 9l (x, Q 2 )p»(S ■ n) + g T (x, Q 2 )S^ + M 2 g 3 (x, Q 2 )n»(S ■ n)] . (1.4) 

From rotational invariance, it follows that the left hand side of ( |1.4|) is proportional to the 
spin vector and thus gi,T,3(x, Q 2 ) must satisfy 

1 rl 

dxgi(x,Q 2 )= dxg T (x,Q 2 ), (1.5) 
•1 J-i 
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jT 1 dx 9l (x,Q 2 ) = 2 dxg 3 (x,Q 2 ). (1.6) 



The same argument for (|1.2j ) leads to the sum rule relations for /ii,z,3(x, Q 2 ): 

'-1 



y dxh\(x,Q 2 )= J dxh L (x,Q 2 ), (1.7) 
y dxh^x^Q 2 ) = 2 y dxh 3 {x,Q 2 ). (1.8) 



The sum rule (|1.5| ) is known as Burkhardt-Cottingham sum rule H and (|1 . 7| ) was first 
derived in Refs. [H, EL H]. Since the twist-4 distributions g 3 , h 3 are unlikely to be measured 



experimentally, the sum rules involving those functions ( Lq) and (|1.8|) are practically useless 
and will not be addressed in the subsequent discussions. Since the left hand side of ( |1.5| ) and 
(|1.7|) are, respectively, the axial charge and the tensor charge, the integral itself in the right 
hand side of ( |1.5| )-( |I~8| ) are finite. As is clear from the above derivation, these sum rules 
are mere consequences of the rotational invariance and there is no doubt in its validity in 
mathematical sense. However, if one try to confirm those sum rules by experiment, a great 
care is required to perform the integral including x = 0. For example, in DIS, x is identified 
as the Bjorken's variable Xb = Q 2 /2P ■ q and x = corresponds to P ■ q — > 00 and this limit 
can never be achieved in the rigorous sense. Accordingly, if hi{x,Q 2 ) has a contribution 
proportional to 5(x) and hi(x,Q 2 ) does not, experimental measurement can never confirm 
the sum rule (|1.7| ) but would rather claim the violation of the sum rule. Such a contribution 
has already been suggested in Refs. [§, |^|. In Refs. || ^ §] it was also argued that the 
integral of Iil(x,Q 2 ) — hi(x,Q 2 ) can be related to the value of certain 'surface terms' that 
appear in formal manipulations involving integrations by parts. 

The purpose of this paper is to reexamine the sum rules involving the first moment of 
the twist-3 distributions. In particular we will argue that the twist-3 distribution Hl{x, Q 2 ) 
has a potential 5-singularity at x = 0, assuming that the twist-2 distributions gi(x,Q 2 ) 
and hi(x,Q 2 ) do not have such singularity. The paper is organized as follows: In section 
2, we examine the sum rule for h^. Starting from the general decomposition of hi based 
on the QCD equation of motion, we will show that it contains a function h™ which has 
<5(x)-singularity (Section 2.1). In Section 2.2, we construct Jil for a massive quark ^from 
the moment of h\ in the one-loop level and show that the h\ also has an 5(x)-singularity, 
which together with the singularity in h™ will give rise to a 5(x) singularity in Hl itself. In 
section 2.3, we will perform an explicit light-cone calculation of hi in the one-loop level to 
confirm the result in the previous sections. Sections 3 and 4 are, respectively, devoted to 
similar examination for the sum rules for e(x, Q 2 ) and gr{x, Q 2 )- 

2 hi(x 1 Q 2 ) for a massive quark 



2 



2.1 <5(x)-functions in hi(x, Q 2 ) 

The OPE analysis of the correlation function ( T71 ) allows us to decompose Hl(x, Q 2 ) into the 
contribution expressed in terms of twist- 2 distributions and the rest which we call h\{x, Q 2 ). 
Since the scale dependence of each distribution is inessential in the following discussion, we 
shall omit it in this subsection for simplicity. Introducing the notation for the moments on 

[-1,1], 



M n \hi 



-i 



dx x n hi J (x), 



this decomposition is given in terms of the moment relation 
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M n [h L ] 

M [h L ] 
with the condition 



n + 2 



M n [h x 



n m r 



n + 2 M 



M n -i[gi 



M n [h 3 L ], (n > 1) 



M [hl] 
Mx[hl] 

By inverting the moment relation, one finds 
h L {x) = hY W (x) + h™(x) + h 3 L (x) 



0. 

: 0. 



2x / dy 



hi(y) 



+ 



m q 
M 



X 



2x I dy 



9i(y) 



+ h 3 L (x) 



2x / dy 



h(y) 



-i 



M 



x 



+ 2x 



(x > 0) 

Mv) 



dy 
i y 

(x < 0) 



+ h 3 T (x) 



(2.1) 

(2.2) 
(2.3) 

(2.4) 
(2.5) 

(2.6) 



(2.7) 



where the first and second terms in Eq. fl2.6| ) denote the corresponding terms in ( |2.7| ). In 
this notation the sum rule fl2.3|) and the condition ( ^4]) implies 



M [h r i 



0. 



(2.{ 



If one naively integrates ( p.7|) over x at x > and x < 0, while dropping all surface terms 
one arrives at Jq dxhi(x) = Jq dxh±(x) + Jq dxh\(x) and likewise for J^ 1 dxhi(x). This 



1 Morc precisely, the original OPE tells us Aia[h^_ 



0. But as long as gi(0±) is finite, which we 



will assume, this is equivalent to stronger relations (2^4) and (|2.8| ) 
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result, in combination with ( |2.4| ), implies that J_ x dxhi = f_ x dxh\. In the following we will 
argue that this procedure may be wrong due to the potentially very singular behavior of the 
functions involved near x = 0. Investigating this issue in detail will be the main purpose of 
this paper. 

We first address the potential singularity at x = in the integral expression for h™(x) 
in ( |2.7| ). In order to regulate the region near x = 0, we first multiply by x", integrate 

from to 1 and let (3 — > 0. This yields 



o+ 



dxKlHx) = ^lim/3 /* dxx^gUy) = ^%i(0+) 



(2.9) 



while multiplying Eq. (|2.7p by \x\^ and integration from —1 to yields 



' — 771 — 771 

i dxhJE{x) = lim/3 d*!^" 1 ^) = -^i(O-), 

(2.10) 

where we have assumed that #i(0+) and gi(0— ) are finite. Adding these results we have 

^ ds^x) + dx^(x) = (^(0+) - ^i(O-)) . (2.11) 

Since we have the sum rule ( |2.8| ) and, in general, lim^o g x (x) — g\(—x) ^ 0, Q we are lead 
to conclude 

777 

W£{x) = K{x) reg - ^± ( gi (o+) - gi ( -)) S (x), (2.12) 

where h™(x) reg stands for the part which is defined by the integral in ( |2.7| ) at x > and 
x < and is regular at x — 0. Since it is unlikely that hY w {x) contains a 5- function at the 
origin, the relation ( |2.12| ) indicates that has a 5(x) term unless h 3 L (x) has a 5(x) term 



and it cancels the above singularity in h™{x). 

Equation ( |2.12| ) clearly demonstrate that the functions constituting hi{x) are more sin- 



gular near x = than previously assumed and that great care needs to be taken when 
replacing integrals over nonzero values of x by integrals that involve the origin. In particu- 
lar, if it turns out that hi{x) itself contains a 5(x) term, then ( |2.3| ) implies 

1 (h L {x) - hix)) + ( (h L (x) - hi(x)) ^ 0, (2.13) 
o+ J-i 

since h\(x) is free from singularity at x = 0: 

0- rl rl 

dxhi(x) + / dxhi(x) = / dxhi(x). (2-14) 
-1 Jo+ J-i 



2 For example, dressing a quark perturbatively at 0(as) yields <?i(0+) ^ and gi(0— ) = gi(0+) = 0. 
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Accordingly an attempt to verify the "/i^-sum rule" || would obviously fail. 

However, as we mentioned earlier, in order to see whether the 5(x) identified in ( |2.12 ) 
eventually survives or not in hi(x), we have to investigate the behavior of h 3 L (x) at x = 0. 
To this end we will explicitly construct hi(x) for a massive quark to 0(a>s)- 



2.2 hi(x,Q 2 ) from the moment relations 

In this subsection we will construct Hl(x, Q 2 ) for a massive quark to 0(a$) from the one-loop 
calculation of A4 n [h 3 L ]. 

One-loop calculation for a massive quark (mass m q ) gives hi(x,Q 2 ) in the following 
form: 



h L (x,Q 2 ) = hf(x) + ^C F \n^h^(x) 



"5. 



Q 2 1.(1), 



2tt 



(2.15) 



where the scale Q 2 is introduced as an ultraviolet cutoff and the Cf = 4/3 is the color factor. 



j^ww,3,m(o,i) ^ defined similarly, ^^(j;) = h^ 1 (x) = 5(1 — x) gives h^'(x) = 5(1 — x), 
as it should. One loop calculation for g\(x) and h\(x) for a quark yields the well known 
splitting functions M, |10[ : 



,(o), 



1 + x 2 3 r/1 

— + -5(1 

[l-x} + 2 1 



X 



^r + 2 5(1 -^ 



(2.16) 
(2.17) 



Inserting these equations into the defining equation in 



one obtains 



hY W{1) (x) 



3x + Axln 



1 — x 



X 



(2.18) 



h 



m(l). 



[1-x) 
2 



1 



x) 



4xln^— ^ - 3 + ^5(1 - x) + ( 3x - ^5(1 



— 4xln- 



x 



3 + 3x 



x 



-5(x) 



-5(x) 



(2.19) 



In the first line of ( 2.19 ), the term (3x — |5(1 — x)) comes from the self-energy correction, 



i.e. from expanding M 

m(l) 



in Eq. O, and -|5(x) = -± gi {0+)6{x) 



in /i™ (x) accounts for the second term on the righthand side of Eq. ( |2.12| ). This term is 
necessary to reproduce the original moment relation ( |2.8| ). We also note that Ji^ W<yl ^ does 



not have any singularity at x = and satisfies Jq 1 dx h L 



WW(1). 



x 



it should. 
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h^\x) can be constructed if we know the purely twist-3 part h^ l \x) in the one-loop 
level. h^(x,u 2 ) can be written in terms of the quark-gluon light-cone correlation function 

as (e = o, r = o) o, n m 



h 3 L (x,Q 2 



ip + 



2tt 



^2ixP+e.- 



udu 



tdt 



x{ps^{-un^i^^gF v a {tnei^{un\ Q APS\\). 



(2.20) 



Starting from this expression, one can, in principle, obtain h^ 1 \x). Alternatively, and 
this is the approach that we will use, one can construct it from already existing one-loop 
calculation of its moments. Taking the n-th moment of (|2.20|) or from the OPE analysis of 
h L M, we have 



doc oc h (^cc . 



[(n+l)/2] 

E 

k=2 



1 - 



2k 



n + 2 2M 



(P5|||i? n , fc (g 2 )|P5||) 



(2.21) 



with 

Rn,k(Q 2 ) 



1 -r 



^(0)a a %z 75 (m ■ D) n ' k igF u a n u (in ■ D) k - 2 ^(0)\ Q 2 - (k -> n - k + 2)(2.22) 



One-loop renormalization of hi was completed in fl4| and the mixing matrix for the local 
operators contributing to the moments of h^x, Q 2 ) was presented. In particular, matrix 
elements of R n> k for a massive quark is given in eq.(3.18) of |L4|]. Since h^\x) for the 
anti-quark is zero, we obtain for the moment of the quark distribution as 



dx x n h 



3(1), 



X 



[(n+l)/2] 

2 E 

1=2 



21 



1 



n 



2U[Z-1] 3 [n-l + 1], 



6 



n + 1 n+2 



(2.23) 



for n > 2 with [fc] 3 = fc(fc + l)(k + 2). [The prefactor in ( [2.23|) is determined by comparison 
with the anomalous dimension for hi and by noting that the operator basis for the quark 
mass operator in [|l4j] has a sign opposite to those for gi] From this result and the defining 
relation for the lowest two moments of h^, (|2.4j ) and fl2.5|) , we can construct (x) as 



hf\x) 



3 - 6x + -8(1 - x) - -5{x) 



(2.24) 



We emphasize that the —l/25(x) in 

3(1) 



is necessary to reproduce the n = moment of 



h% '(x). From ( pA8|) , ( pJ9|) and ( ggg ), one obtains 



^L(x,Q 2 ) = 5(l-a;) + ^ln^C F 

27r 



[1 



^p + ^Cl-x)-^) 



(2.25) 



We remark that the above calculation indicates that the 8{x) term appears not only in h™ 
but also in h\. Furthermore they do not cancel but add up to give rise to —S(x) in hi(x, Q 2 ) 
itself. 

In the next subsection we will confirm Eq. (|2.25| ) through a direct calculation of hi(x, Q 2 ) 
for a quark. 



2.3 Light-cone calculation of hi{x,Q 2 ) 

In order to illustrate the physical origin of the 8(x) terms in h^x), and in order to develop 
a more general and convenient procedure for calculating such terms, we shall now evaluate 
Iil(x) using time-ordered light-front (LF) perturbation theory. The general method has 
been outlined in Ref. |16| and we will restrict ourselves here to the essential steps only. 

There are two equivalent ways to perform time-ordered LF perturbation theory: one 
can either work with the LF Hamiltonian for QCD and perform old-fashioned perturbation 
theory (the method employed in Ref. JRJ), or one can start from Feynman perturbation 
theory and integrate over the LF-energy k~ first. We found the second approach more 
convenient for the one- loop calculation of Hl{x) and this is what we will use in the following. 

In LF gauge, A + = 0, parton distributions can be expressed in terms of LF momentum 
densities (A; + -densities). Therefore, one finds for a parton distribution, characterized by the 
Dirac matrix T at 0(as) and for < x < 1 



fr(x)u(p)Tu(p) = -ig u(p) 



d A k ( k + ' 



(2tt 



r 



— m q + is 



m„ + is 



l v u(p)D^{p- k), 



where 



(2.26) 



DM 



1 



q l + is 



9nv 



q^n u + n^q u 
qn 



(2.27) 



is the gauge field propagator in LF gauge, and n M is a light-like vector such that nA = A + ~ 
(A + A 3 ) /y/2 for any four vector A>*. 

The k~ integrals in expressions like Eq. ( |2.26| ) are performed using Cauchy's theorem, 
yielding for < k + < p + 



— i 



dk 



1 



2tt 



(k 2 



m 2 + is 



y(p-k) 2 + ie (2k+) 2 2(p+ - k+) ^_ 



2k+ 2(p+-k+) 



1 1 



X 



2^ 



(2.28) 



where we used k + = xp + . In order to integrate all terms in Eq. ( |2.26| ) over k , Cauchy's 
theorem is used to replace any factors of k~ in the numerator of Eq. (|2.26|) containing k~ 
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by their on-shell value at the pole of the gluon propagator 



k = p 



IP-L 



(2.29) 



In the following we will focus on the UV divergent contributions to the parton distribution 
only. This helps to keep the necessary algebra at a reasonable level because we restrict 
ourselves to the leading behavior in k^, which arises from those terms in the numerator of 
Eq. (|2.26|) that are quadratic in k^, and therefore give rise to a logarithmically divergent k^ 
integral. The transverse momentum cutoff can be replaced by Q 2 in these expressions. The 
rest of the calculation is just tedious algebra and we omit the intermediate steps here. 
We find for < x < 1 to 0(a s ) 



h r (x,Q 2 ) = ^C F \n Q 



L \x, Q ) - —C F „„ 

27r m A q [1 — x\ 



(2.30) 



1, i.e. 



-J— for x < 1 and 

l—x 



where the usual +-prescription for n^i applies at x ...... ^ 

f dx ^\ = 0. Furthermore, hz(x) = for x < 0, since anti-quarks do not occur in the 
0(ag) dressing of a quark. In addition to Eq. (|2.30| ), there is also an explicit 6(x — 1) 
contribution at x — 1. These are familiar from twist-2 distributions, where they reflect the 
fact that the probability to find the quark as a bare quark is less than one due to the dressing 
with gluons. For higher-twist distributions, the wave function renormalization contributes 
is f^-CV In — 1). The same wave function renormalization also contributes at twist-3. 

However, for all higher twist distributions there is an additional source for 8{x — 1) terms 
which has, in parton language, more the appearance of a vertex correction, but which arises 
in fact from the gauge-piece of self-energies connected to the vertex by an 'instantaneous 
fermion propagator' For gx(x,Q 2 ) these have been calculated in Ref. fl6|] where they 

give an additional contribution — ffrCp In ^S(x — 1), i.e. the total contribution at x — 1 

for gr(x, Q 2 ) was found to be f^-Cp In ^\5{x — 1). We found the same 5(x — 1) terms also 
for h^x, Q 2 ). [] Combining the 5(x — 1) piece with Eq. ( |2.30| ) we thus find 



h L (x,Q 2 ) = 5{x -1) + ^C F ln^ 

Z7T Tflq 



r 2 : +-5(X-1) 

[l-x} + 2 1 J 



for x > 0. (2.31) 



Comparing this result with the well known result for hi [110 
hl ^Q^ = S(x-l) + ^\n%C F 2r 



2tt 



[1 



x\ 



+ -5(x-l) 



(2.32) 



3 Physically, the reason why the wave function renormalization contribution depends on the twist is 
that in LF gauge, different components of the fermion field aquire different wave function renormalization. 
However, since all twist-3 parton distributions involve one LC-good and one LC-bad component, it is natural 
to find the same wave function renormalization for all three twist-3 distributions. 
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one realizes that 



lim f 1 dx \h L (x, Q 2 ) - h(x, Q 2 } = — In %C F + 0, (2.33) 
£^oj £ l J 27T m z q 

i.e. if one excludes the possibly problematic region x — 0, then the /i^-sum rule pi is 
violated already for a quark dressed with gluons at order 0{as). 

In the above calculation, we carefully avoided the point x = 0. For most values of k + , 
the denominator in Eq. Q2.26] ) contains three powers of k~ when k~ — > oo. However, when 
/c + = 0, A; 2 — m 2 becomes independent of and the denominator in Eq. (|2.26|) contains 
only one power of k~ . Therefore, for those terms in the numerator which are linear in k~, 
f] the ^"-integral diverges linearly. Although this happens only for a point of measure zero 
(namely at k + = 0), a linear divergence is indicative of a singularity of h^ix, Q 2 ) at that 
point. To investigate the k + pa singularity in these terms further, we consider 

f(k + ,k ± ) EE J dk~ 



dk~ 



(k 2 — m 2 q + ie) 2 {p — k) 2 + ie 
k~ + ( k~ - V 



(k 2 — m 2 + ie) 2 [(p — A;) 2 + ie] 

= / can .(^ + ,k±) + /^n.(^ + ,k ± ), (2.34) 

where the 'canonical' piece f can . is obtained by substituting for k~ its on energy-shell value 
kr (f2~29|) (the value at the pole at (p - k) 2 = 0) 

For k + = xp + 7^ 0, it is only this canonical piece which contributes. To see this, we note 
that 

k--k- = - n j P ~ k }' (2.36) 



and therefore 

/ sin (fc + ,k ± ) = J dk' 



2(p+ -k+) 



kr - kr 



{k 2 — m 2 + ie) 2 (p — k) 2 + ie 



2(p+-k+)J dk (k 2 -m 2 + ie) 2 ' (2 ' 37) 



4 This is the highest power of k that can appear in the numerator for twist three distributions. 
5 Therc is another point, k + = p + , where a similar divergence occurs, but the latter is only logarithmic. 
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Obviously |TT| 



/ 



dk 



(2k+k- 



k 2 ± — m 2 + is) 



(2.38) 



m„+k| —is 



for fc + 7^ because then one can always avoid enclosing the pole at k~ = " Lq ^i by 
closing the contour in the appropriate half-plane of the complex k~ — plane. However, on 
the other hand 



dk dk 
and therefore 



(2k+k~ - k\ - m 2 + isf 



d 2 k, 



1 



Z7T 



ik\ — — m 2 + is) 2 k\ + m\ 



fsin{k 7 kj 



1 m5{k^ 



2p~ 



k 2 , +m 2 ' 



(2.39) 



(2.40) 



Upon collecting all terms oc k in the numerator of Eq. Q2.26|) , and applying Eq. (|2.40|) to 
those terms we find after some algebra those terms in hi{x, Q 2 ) that are singular at x = 



h L ,sin(x,Q 2 ) = - < ^\n^-C F 5(x). 

Z7T m A 



(2.41) 



Together with Eq. (|2.31|) , this gives our final result for h L , up to 0(a s ), valid also for x = 



h L (x, Q 2 ) = 5(x 



1 + — C^f hi - 

2tt m l q 



-5{x) + 



[1-x] 



+ -S(x-1) 



(2.42) 



As expected, Hl from Eq. ( |2.42| ) does now satisfy the /i^-sum rule, provided of course the 
origin is included in the integration. 

This result is important for several reasons. First of all it confirms our result for hi(x, Q 2 ) 
as determined from the moment relations. Secondly, it provides us with a method for 
calculating these 6(x) terms and thus enabling us to address the issue of validity of the 
naive sum rules more systematically. And finally, it shows that there is a close relationship 
between these 5(x) terms and the infamous zero- modes in LF field theory [fLgfl . 

While we were completing the manuscript for this paper, we learned of Ref . |15| , where 
canonical Hamiltonian light-cone perturbation theory is used to calculate hi(x). For x ^ 
the result obtained in Ref. |T5| agrees with ours which provides an independent check of the 



formalism and the algebra. However, the canonical light-cone perturbation theory used in 
13 1 is not adequate for studying the point x = 0. From the smooth behaviour of fiL{x) 
the authors of Ref. |15j conclude that the sum rule for the parton distribution 



Ref. 



near x 

Hl{x) is violated to 0(as)- Our explicit calculation for h^x) not only proves that the sum 
rule for Hl{x) is not violated to this order if the point x = is properly included, but also 
shows that it is incorrect to draw conclusions from smooth behaviout near x = about the 
behaviour at x — 0. 
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3 e{x,Q 2 ) 



The other chiral-odd twist-3 distribution e(x, Q 2 ) is also expected to satisfy a simple operator 
sum rule 

1 dxe(x,Q 2 ) = J-(P\ VS(0)V(0) ^ \P), (3.1) 



i 



2M 



which follows trivially by integrating (|1.3|) over x. Because of our results from above for 
Hl{x, Q 2 ), we are now of course more cautious and address in the following the issue whether 
Eq. (|3.1|) is also valid if the origin is excluded from the region of integration. Again, we will 
consider a massive quark to 0(as) 

3.1 Constructing e(x, Q 2 ) from its moments 



The OPE analysis of (|1.3| ) decomposes e(x) into the twist-2 contribution and the purely 
twist-3 piece e 3 (x) as 

777 

M n [e] = M n [e 3 ] + j^M n -i[fi], (n > 1) (3.2) 
with the relation (|3.1| ). This moment relation defines the decomposition 

777 

xe(x, Q 2 ) = xe 3 (x, Q 2 ) + ^A(x, Q 2 ). (3.3) 

where we multiplied x to regularize any possible 6(x) contributions, which will be specified 
later. The one loop calculation of e(x, Q 2 ) for a massive quark yields 

e(x, Q 2 ) = 5(1 -x) + ^C F \n^-e^(x). (3.4) 

27r 

The lowest moment of (x) can be obtained directly from the cr-term relation ( ^.1[ ) 

MM = = = ™ = ! + %C F x (3.5, 



i.e. the lowest moment of e'- 1 '' reads 



Corresponding to (|3.2j) , we have 



M [e^} = | (3.6) 



M n [e®] = M n [eW) + [K/i) (1) ] , (n > 1) (3.7) 
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Our problem is to construct e^'(x) for a massive quark from ( ^.6| ) and ( |3.7| ). In principle 
we can obtain the purely twist-3 piece e 3 ^ from the one-loop calculation starting from the 



correlation function [n], 0, [13 



J2ixP-i 



2tt 



d« / rft(P|^(-<)^ Q ^ Q (^)U^«)lQ^)- (3-8) 



However, we again make a short cut to get it from the moment. The moment of e 3 W(x) 
for a massive quark is given in (19| as a part of the mixing matrix in the context of the 
renormalization. From eqs.(4.2) and (3.18) of [19] we have for the n-th moment of e 3 ^ 
(n>l): 



-i 



dxx n e 3{1 \x) 



n/2 
1=2 L 



[I - 1] 3 [n-l + 1], 



[n/2], 



(n+l)/2 

2 E 

Z=2 



+ 



[I - 1] 3 [n-l + 1] 3 



1 1 1 

« + 7 

2 n n + 1 



which gives 



(n : even) 
(n : odd) 



(3.9) 



xe 3 ^(:r) = x-5(l — x) — 1 + x. 



(3.10) 



Together with the twist-2 contribution 



1 

M 



™ji) (1) 



one obtains 



1 + x 2 

[i - 4 



[l-x]+ 2 1 J 



for x > 0. 



(3.11) 



(3.12) 



Note that 1/x singularity in e 3 ^(x) and /] /x cancel each other in and e^(a;) itself 
is integrable at x = 0. In order to satisfy the cx-term sum rule ( |3.6|) one needs to introduce 
5(x) term. Accordingly the final result, which satisfies the moment relations both for n = 
and n^O, reads 



(3.13) 



This result clearly indicates that the sum rule (3.1) is satisfied only by including the point 
x = in the integration. 
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3.2 Light-cone calculation 

We start from Eq. ( ggg ) with T = 1 

d 4 k 



e(x)u(p)Tu(p) = —ig u(p) 



5 x 



r 



— m q + ie tp — m q + ie 



l v u{p)D^ u {p-k). 

(3.14) 



the numerator algebra in Eq. ( 3.14| ) yields 



/IV 



(p - k) u n^ + (p- kyn v 



P 



2u{p) 



k +m 2 + 2 



k 2 p + — m 2 k + 



2m n 



k+ 



u{p). 



(3.15) 



p + — k + 

In order to determine the 'canonical' part of e(x, Q 2 ) we use again contour integration, 

picking up the pole at k~ = k~ for < k + < p + , which allows us to replace k 
k 2 



2 . ill 
1-x 



andp+A; -> 2(1 _ x) , 
yields for < x < 1 



where we kept only the leading terms in kj_. After some algebra this 



e w (x) 



1 — x 



(3.16) 



in agreement with the result obtained from the moments. For the wave function renormal- 
ization (the coefficient in front of 5(x) the same coefficient is obtained as for Hl(x, Q 2 ) and 
the details will be omitted here. 

Finally, we focus on possibly singular terms near x = 0. The only numerator term 
involving k~ in Eq. ( 3.15Q , which is not multiplied by k + appears in — 4m q u(p)k'u(p). Upon 
repeating the same steps as in Section 3, i.e. isolating the singular piece by adding and 
subtracting k~ , canceling k~ — k~ against the gluon propagator, and making use of Eq. 
fl£30D one finds 



Collecting all terms we thus find 



e (1) fa;) 



e {1 \x)sin = S(x). 



^+^(1 -*)+*(*)., 



(3.17) 



(3.18) 



which completely agrees with the result from above Q3.13j ) 
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4 g T (x, Q T 
4.1 Possibility of 5(x) in gx 

From a practical point of view, gr{x, Q 2 ) is the most important among the twist-3 distribu- 
tions, because it is the least difficult to measure experimentally. For this reason, gx^x^Q 2 ) 
has been the subject of many studies in the literature. 

We again start with the decomposition of gr(x) into the Wandzura-Wilcek term g^ w (x) 
20], the term proportional to the quark mass and the purely twist-3 part g^(x): f\ 



9t{x) 



g? W (x)+g?(x)+gUx) 



y M 



hi(x) 



x 



X dy 9 -^ 

-i y 



m q 
M 



hi(x) 



+ 



x 



-i 



dy 



hi(y) 



(x > 0), 



(x < 0). 



(4.1) 



(4.2) 



This decomposition is obtained by the analysis of the correlation function using QCD equa- 
tion of motions. The corresponding moment relations are 



M n [g T ] 
M [g T ] 



1 



n + l 
M \gi], 



M n -i\gi 



n m q 
n+l M 



M n -![hi] + M n [g%], 



(n > 1) 



with 



M [g 3 T }=0, 



where we assumed that g^(x) itself is integrable at x = 0. The relation 
the Burkhardt-Cottingham (BC) sum rule0. We note ( |4.4| ) and fl4.5| ) implies 



(4.3) 
(4.4) 



(4.5) 
is known as 



M \g\ 



0. 



(4.6) 



In the above discussion we again assumed separate relations ( [4.5D and ft4.6|) instead of 
Mo[gT + = 0. This is justified as long as /ii(0±) is finite. (See below.) 
Following the same argument leading to (|2.12|) for h™(x), we have 



This relation together with 

g ™{x) =gf[x)\ reg 



dxgjp(x) + / dxg?{x) = -j (/n(0+) - ^(0-)) . 

implies g™{ x ) has a singularity at x = as 



M 



(/n(o+)-/n(o-))<y(x) 



(4.7) 



(4- 



6 Here and below in this subsection we omit the scale dependence of the distribution functions for 
simplicity. 
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where g™(x)\ reg is a part obtained from the integral ( |4.2p at x > and x < and is regular 
at x = 0. This relation shows that if hi(0—) = — /ii(0+) 7^ /ii(0+) then the Burkhardt- 
Cottingham relation would be violated if data is taken only for nonzero x — unless of course 
there is another 5(x) contribution to g^(x) which happens to cancel exactly the one in g™. 
Perturbation theory predicts at small x that h\ (x) ~ x in the leading order as is seen from 
( |2.17D and hi(x) — hi(—x) ~ const, at the next-to-leading order [2T] . BFKL approach gives 



h\(x) ~ const, at small x [p2|| . Therefore there is a possibility that the 5{x) term in 
survives and hence a seeming violation of BC sum rule. 

In the next subsections, we will calculate gr(x, Q 2 ) for a massive quark to 0(as) in order 
to look for (possibly) another origin of the 5(x) contribution. 

4.2 Constructing gr(x, Q 2 ) from its moments to 0{cts) 
One loop calculation of gr{x) for a massive quark takes the form of 

g T (x, Q 2 ) = 6(1 -x) + ^C F \n^g?(x). (4.9) 

Z7T mi 



Corresponding to the decomposition ( |4.1|) we write 

g?(x)=g^ 1 \x) + gp 1 \x) + g^\x). (4.10) 

Inserting the expression ( |2.16| ) for g± and ( |2.17| ) for into ( [4.2| ), and also taking into 
account the self-energy correction, we have 

g™ W{1 \x) + gp\x) = \nx+ 1 - + x + (4.11) 

1 [i — x\ + 

This result has no 5(x) term, which is simply because hi(0 + ) = hi(0~) = in the one-loop 
calculation. (See ([4.81).) 

As in the case of h\(x) and e 3 (x), we construct from its moment for a massive 

0/-1 \ 

quark. From eqs.(5) and (18) of ||23|| , we have for the n-th moment of g T ; (x) as 

1 9. 1 

+ T-r^- ( 4 - 12 ) 



2 2(n+l) (n + 1) 2 ' 

This moment relation is originally given for n > 1, but it also gives zero for n — 0, which is 
consistent with ( f4.5| ). Accordingly no S(x) term is necessary to obtain 

gf\x) = U{l-x)-^-\nx. (4.13) 



From (|4.11| ) and (|4.13|) , we get 



g T ( x ,Qi) = 5(l-x) + ^C F \n®- 



l + 2x-x 2 1 

H (3(1 — xj 



2tt [ [1 - x}+ 2 



(4.14) 
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4.3 LC calculation 



We calculated gr(x, Q 2 ), using using the same LF pQCD techniques that we used to study 
Hl{x,Q 2 ) and e(x, Q 2 ). The algebraic steps involved are rather lengthy and we therefore 
present only the final result here, which reads 



g T (x,Q 2 ) = 6(l 



x) + — Cp In — - 



2tt 



mi 



1+x 2 l r/ 



(4.15) 



161, which 



without any 5(x) term. This result completely agrees with the findings from Ref. 
confirms our formalism. 

Even though the numerator for T = 7^75 in Eq. (|2.26|) contains k~, those terms turn 
out to be multiplied by at least one power of k + . Since x5(x) = 0, there are no 5(x) terms 
in gx{x,Q 2 ) to 0(as)- However, we do not have a simple explanation (other than working 
out the numerator algebra) as to why factors of k~ in the numerator algebra are always 
accompanied by at least one power of k + for gx^x.Q 2 ) and not for hi(x,Q 2 ) and e(x, Q 2 ) 
in this one loop calculation. 



5 Summary 



We have investigated sum-rules for twist-3 distributions in QCD, and found examples where 
these sum-rules are violated if the point x = is not properly included. 
For a massive quark, to 0(as) we found 



g T (x,Q 2 ) = S(x 
h L (x, Q 2 ) = 5(x 
e(x, Q 2 ) = S(x 



D + 



D + 





In 












In 




2tt 




m 2 


^C F 


In 




2tt 




m 2 



1 + X * 1 r/ 

2x+ [T^xU + 2 Six 



-S(x) + 



S(x) + 



[1 



[1 



X 



+ ^(x - 1) 



+ -5(x-l) 



(5.1) 



At 0(as) neither hi(x, Q 2 ) nor e(x, Q 2 ) satisfy their respective sum rule if one excludes the 
origin from the region of integration (which normally happens in experimental attempts to 
verify a sum rule). gT{x,Q 2 ) is the only exception and its sum- rule is satisfied even when 
the origin is not included. 

Of course, QCD is a strongly interacting theory and parton distribution functions in 
QCD are nonperturbative observables. Nevertheless, if one can show that a sum rule fails 
already in perturbation theory, then this is usually a very strong indication that the sum 
rule also fails nonperturbatively (while the converse is often not the case!). 



16 



From the QCD equations of motion, we were able to show nonperturbatively that [] the 
difference between hi(x, Q 2 ) and h 3 L (x, Q 2 ) contains a 5 function at x = 

h L {x,Q 2 )-h 3 L (x,Q 2 ) 

Since gi(0+,Q 2 ) — gi(0—,Q 2 ) = lim^o gi(x, Q 2 ) — gi(x.Q 2 ) seems to be nonzero (it may 
even diverge^), one can thus conclude that either Iil(x, Q 2 ) or h 3 L (x, Q 2 ) or both do contain 
such a singular term. 

We checked the validity of this relation to 0(as) and found that, to this order, both h\ 
and h,L contain a term oc 5(x). We also verified that even though the sum rule for h^x) 
and e(x) are violated if x = is not included, the sum rules for all three twist 3 parton 
distributions are still satisfied to 0(ots) if the contribution from x = (the 5(x) term) is 
included. 



sin.mi.ln.r 9 A/7 



(gi(0+,Q 2 )-gi(0-,Q 2 ))5(x). (5.2) 
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